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ABSTRACT 



Aims. We study the efficiency of Joule heating in the crustal layers of young neutron stars. 

Methods. The induction equation is solved numerically for a realistic crust model with a conductivity that depends on the density 
and temperature. 

Results. We show that dissipation of the magnetic field is highly inhomogeneous in the crust with much faster dissipation in relatively 
low density layers. In young neutron stars, the rate of Joule heating in the crust can exceed the standard luminosity of non-magnetic 
stars and can even be comparable to the luminosity of magnetars. 

Conclusions. The results of calculations are compared with the available observational data. We argued that the crustal field model 
can well account for the data on the surface temperature and magnetic field of young neutron stars. 
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The influence of the magnetic field on the thermal history of neu- 
tron stars is twofold. On one hand, the field changes the character 
of heat transport providing some anisotropy to kinetic processes. 
Such anisotropic transport can result in a characteristic temper- 
ature difference between the poles and equator caused by large- 
scale inhomogeneities of the magnetic field (see, e.g., Potekhin 
et al. 2003; Perez-Azorin et al. 2006). On the other hand, dissipa- 
tion of the magnetic field provides an additional heating mecha- 
nism that can substantially alter the cooling evolution (Miralles 
et al. 1998). Even more or less standard for radiopulsars mag- 
netic fields ~ 10^^ — 10^'^ G can maintain the surface tempera- 
ture as high as ~ 10^ K for ~ 10 — 100 Myr which is substan- 
tially higher than predicted by various cooling models. A gen- 
eral relativistic treatment of Joule heating in neutron stars un- 
dertaken by Page et al. (2000) led to results very similar to those 
of Miralles et al. (1998) since, as it was argued by Konar (2002), 
general relativistic corrections do not significantly influence the 
magnetic field decay in neutron stars. 

The internal magnetic fields in neutron stars, however, can 
be substantially stronger than the "standard" surface field of 
radiopulsars, and they can produce much more efficient heat- 
ing. The nature of pulsar magnetism is still the subject of de- 
bate. One of the possibilities is associated with the turbulent dy- 
namo action that can amplify the magnetic field during the first 
~ 30 — 40 s of a neutron star's life when the star is subject hy- 
drodynamic instabilities (Thompson & Duncan 1993; Bonanno 
et al. 2003). Turbulent motion caused by instabilities in com- 
bination with rapid rotation that seems to be almost inevitable 
in protoneutron stars make a turbulent dynamo one of the most 
plausible mechanisms of pulsar magnetism. The characteristic 
feature of a turbulent dynamo in stars is that generation of the 
large-scale field is always accompanied by generation of small- 
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scale fields that are approximately in equipartition with the tur- 
bulent motion. Therefore, if the field of pulsars is of dynamo 
origin, it should have a very complex geometry when dynamo 
action stops (see Urpin & Gil 2004, Bonanno et al. 2005, 2006). 
Magnetic fields of various scales will be frozen into the crust that 
begins to form soon after the end of the unstable stage and then 
will decay due to the finite crustal resistivity. 

Miralles et al.(1998) have considered the effect of Joule heat- 
ing in old neutron stars with an age t > 10^ yrs, when the surface 
temperature is low and dissipation of even standard pulsar mag- 
netic fields can change the thermal evolution substantially. In this 
paper, we treat the effect of dissipation during the early evolu- 
tion of neutron stars (t < 10^ yrs), when neutron stars are much 
hotter The crustal conductivity depends on the temperature and 
can be relatively low in hot neutron stars. Due to this, the field 
decay is rapid in young neutron stars with an age ^ 10** — 10*^ 
yrs. This type of behaviour has been well known since the paper 
by Urpin & Muslimov (1992). The fast decay has been studied 
in detail by a number of authors (see, e.g., Urpin & Konenkov 
1997; Sengupta 1998; Mitra et al. 1999, etc.), and was rediscov- 
ered with no citation in recent calculations by Pons & Geppeit 
(2007). However, even such efficient field decay during the ini- 
tial stage can influence the thermal evolution only if the magnetic 
field is sufficiently strong and, therefore, our study primarily ad- 
dresses high magnetic field pulsars. 



2. The internal magnetic fields in pulsars 

The unstable stage lasts only ~ 30 — 40 s in proto-neutron stars, 
but this is sufficient for dynamo action to generate a strong mag- 
netic field because the turbulent motion is very fast. A turbulent 
dynamo can generate a large-scale magnetic field only if the star 
rotates rapidly and the rotation period Pq is shorter than some 
"critical" value Per (^ 1 s in the simplest model). Most proto- 
neutron stars likely rotate faster and, as a result, a dynamo should 



2 



V.Urpin and D.Konenkov: Joule heating in high magnetic field pulsars 



operate in these objects (Bonanno et al. 2003, 2005). The dy- 
namo generates both the toroidal field that is essentially internal 
and vanishes at the surface and the poloidal field that is non- 
vanishing outside the star. The poloidal field is responsible for 
pulsar activity. Typically, the generated toroidal field is compa- 
rable to the poloidal one if the proto-neutron star rotates rigidly, 
but it can be substantially stronger if rotation is differential. For 
example, the toroidal field is ~ 20 — 60 times stronger than the 
poloidal one if rotation at the center is three times faster than 
that at the surface equator. The saturation field generated by the 
dynamo can be estimated by making use of the simplest model 
of a-quenching (Bonanno et al. 2006). This yields for the gener- 
ated large-scale field 

Bo ~ B,,g^P„./Pn - 1, (1) 

where Bf.q is the equipartition small-scale field, B^q « Airpv^ 
where vt is the turbulent velocity. The equipartition field varies 
during the unstable stage, rising rapidly soon after collapse and 
then decreasing when the temperature and lepton gradients are 
smoothed. At the peak, B^q ^ 10^® G in the convective zone 
(Thompson & Duncan 1993) and B^g ~ (1 - 3) x 10^^ G in 
the neutron-finger unstable zone (Urpin & Gil 2004). However, 
Beq decreases and becomes smaller than the peak value as the 
protoneutron star cools down. The dynamo operates while the 
quasi-steady condition TcooI S> t is fulfilled with TcooI being 
the cooling timescale and r the turnover time of turbulence. The 
final strength of a generated small-scale magnetic field is ap- 
proximately equal to B^q at the moment when the quasi-steady 
condition breaks down and r ~ TcooI- If t ~ Tcooh we have 
Vt ~ t^^t/tcooI where (.t is the turbulent lengthscale. Then, 

Beq ^ \/ ^TTfWT ~ TT AtT ptxT^oll- (2) 

Estimate (2) yields Beq ^ (1 — 3) x 10^^ G for the maximum 
possible lengthscale It ^ ^ — 5 km and TcooI ^ few seconds. 
Substituting estimate (2) into expression (1), we obtain 

Bo « .[^1 ~ (1-3) X 10-./^ G.(3) 

Tcooi V y Fq 

The shortest possible period is likely ^ 1 ms and, hence, the 
strongest internal magnetic field generated by a dynamo can 
reach ~ (0.3 - 1) x 10^^ G. The field at the surface is usually 
weaker than the internal field generated by a dynamo (Bonanno 
et al. 2006). Note that generation is more efficient in the outer 
layers because the parameter a that determines the efficiency of 
the dynamo is proportional to O and the density gradient (see, 
e.g., Rtidiger & Kitchatinov 1993) but the latter is maximal in 
the outer layers. 

If Pq is longer than Per, then the large-scale dynamo does 
not operate in proto-neutron stars but the small-scale dynamo is 
still efficient. One can expect that such neutron stars have only 
small-scale fields with the strength B^q ^ (1 - 3) x 10^^ G 
and have no dipole field. Likely, slow rotation with Pq > Per is 
rather difficult to achieve in proto-neutron stars, and the number 
of such exotic objects is small. 

It is difficult to predict the final disposition of the mag- 
netic field after turbulent motion stops. The magnetic configu- 
ration generated by the dynamo can be unstable, for example, if 
it contains a strong toroidal component. MHD instabilities can 
lead to a rapid evolution on the Alfven time scale and form a 
configuration with comparable poloidal and toroidal fields (see, 
e.g., Braithwaite & Spruit 2006; Braithwaite & Nordlund 2006). 
Therefore, it is unlikely that the toroidal field is substantially 



stronger than the poloidal one in neutron stars. When convec- 
tion is exhausted, both large- and small-scale fields evolve under 
the influence of ohmic dissipation and buoyancy. Likely, for- 
mation of the crust provides the most important influence on 
the evolution of magnetic fields at this stage. Soon after con- 
vection stops, the neutron star cools down to the internal tem- 
perature ~ (1 — 3) X 10^° K at which neutrons and protons 
can form nuclei and clusters in matter with a density ~ 10^^ 
g/cm^. When the temperature decreases progressively, nuclei 
can be formed at a lower density as well. Nuclear composi- 
tion at high density is subject to many uncertainties and de- 
pends generally on the pre-history of the neutron star. For in- 
stance, the composition of the "ground state" matter (Negele & 
Vautherin 1973) differs noticeably from that of the "accreted" 
matter (Haensel & Zdunik 1990). The Coulomb interactions of 
nuclei leads to crystallization that occurs when the ion coupling 
parameter r = Z'^e^/afc^T reaches the critical value Fm w 170 
(Slattery et al. 1980); a = (3/47rni)^/^ is the mean inter-ion 
distance, rij and Z are the number density and charge number of 
ions, respectively; T is the temperature, and fcs is the Boltzmann 
constant. Then, the crystallization temperature is 

Trn « 1.3 X lO^Z-'^/'^/i.-i/'V}^' K, (4) 

where pe is the number of baryons per electron, and pi2 = 
p/10^^ g/cm^. For p ~ 10^^ - 10^^ g/cm^, the crystallization 
temperature is ~ 10^" K. Therefore, the crust formation starts al- 
most immediately after the end of the convective phase (or even 
before the end), and the magnetic field generated by the dynamo 
should be frozen into the crust. 

3. Joule heating in tlie crust 

Likely, the evolution of the magnetic field in the crust of hot 
neutron stars is mainly determined by the ohmic dissipation, 
and Joule heating has an important influence on the thermal 
evolution. The rate of Joule heating depends both on the ge- 
ometry of the magnetic field and the conductive properties of 
the crust. In a strong magnetic field, the ohmic dissipation can 
be accompanied by Hall currents (Shalybkov & Urpin 1997, 
Hollerbach & Riidiger 2002, 2004). These currents are non- 
dissipative since they do not contribute directly to the rate of 
Joule heating. However, the Hall currents couple different modes 
and alter the magnetic configuration by redistributing the energy 
among modes. In this way, the Hall currents can indirectly affecr 
the rate of dissipation. Note also that, because of the non-Unear 
nature of the Hall currents, they can generate magnetic features 
at the stellar surface with a smaller lengthscale than the back- 
ground magnetic field. This has been known since the paper by 
Naito & Kojima (1994) (see also Muslimov 1994, Shalybkov & 
Urpin 1997) who argued that some neutron stars may have strong 
disordered magnetic fields at their surface due to the HaU effect. 
The same conclusion recently was reached by Pons & Geppert 
(2007). 

The influence of the HaU effect on the magnetic field is char- 
acterized by the HaU parameter, w^r, where lob = eB/m^c 
is the gyrofrequency of electrons and r is their relaxation time, 

is the electron effective mass. The Hall effect is unimpor- 
tant if ujbt < 1, but it can influence the magnetic evolution if 
the field is sufficiently strong and ojbt ^ 1. Indeed, numerical 
simulations (Shalybkov & Urpin 1997, Hollerbach & Riidiger 
2002, 2004) indicate that some acceleration of the decay can oc- 
cur, particularly if the Hall parameter is very large, ojbt > 200. 
This conclusion, however, turns out to be sensitive to the initial 
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magnetic geometry. For a purely toroidal field, the Hall effect 
develops fine field structures, which cause the field to decay on 
the Hall time-scale rather than on a much longer ohmic time- 
scale. However, the presence of even a relatively weak poloidal 
field (~ few per cent) changes the result drastically, and dissipa- 
tion proceeds approximately on the ohmic timescale. Since the 
magnetic field in the neutron star should have a rather complex 
geometry with both the poloidal and toroidal components, we 
can expect that departures in the rate of Joule heating caused 
by the Hall effect have to be small compared to the case when 
the ohmic dissipation is acting alone. Note that some simplified 
models describe the influence of the Hall effect on the magnetic 
evolution in terms of the decay on the Hall timescale (see, e.g.. 
Gumming et al. 2004). From a comparison with numerical mod- 
elhng, it is seen that this naive picture is insufficient. The Hall 
timescale often characterizes the period of oscillations rather 
than the decay timescale of magnetic modes. Such an oscillatory 
behaviour of the decaying magnetic modes was first obtained by 
Urpin & Shalybkov (1997). Despite Pons & Geppert (2007) hav- 
ing confirmed the earher finding by Shalybkov & Urpin (1997) 
and HoUerbach & Riidiger (2002, 2004) tiiat the HaU effect does 
not substantially accelerate dissipation of realistic magnetic con- 
figurations, Aguilera, Pons, and Miralles (2007) in their study of 
Joule heating in neutron stars assume confusingly that the Hall 
currents result in a fast decay of the magnetic field. 

The field decay can also be influenced by the so-called Hall- 
drift induced instabiUty if it occurs in neutron stars (Geppert & 
Rheinhardt 2002). These autiiors found tiiat instabiUty can occur 
which raises small-scale magnetic modes if the initial magnetic 
field is strong and its second spatial derivative is large enough. It 
is unclear, however, whether the effect of this instabiUty on the 
dissipation rate is noticeable or not since a non-linear regime has 
not been considered. In particular, HoUerbach & Riidiger (2002) 
express doubts that a linear instability starting from small am- 
plitude perturbations can be compatible with the cascade caused 
by the non-linear Hall effect. Since even this cascade leads only 
to smaU departures in the dissipation rate for realistic magnetic 
configurations (Shalybkov & Urpin 1997, HoUerbach & Riidiger 
2002, 2004), the effect of a linear instability should be very 
weak, if it exists. As far as the non-linear regime of instabil- 
ity is concerned, Geppert et al. (2003) assume that, at saturation, 
the ampUtude of unstable perturbations can reach a significant 
fraction of the background field and, hence, the latter should de- 
cay on the Hall timescale. Since this obviously contradicts the 
observed lifetime of pulsars (~ 10^ yrs), the authors assume 
that the Hall drift induced instabUity can operate only during 
part of the neutron star Ufe whereas decay before and after such 
episodes is determined by the ohmic dissipation. According to 
Rheinhardt & Geppert (2002), episodes occur when the age of a 
neutron star is > 10^ yrs since the Hall parameter is not suffi- 
ciently large in younger stars. However, studies of pulsars based 
on methods of population synthesis (Bhattacharya et al. 1992, 
Hartman et al. 1997, Regimbau & de Freitas Pacheco 2001, 
Faucher-Giguere & Kaspi 2006) indicate that the magnetic field 
of radiopulsars decays Uttle during their active lifetime. A signif- 
icant magnetic field decay in radiopulsars Ukely can occur only 
in very young stars with age < 10^ yrs since the number of these 
pulsars is small and they give a negligible contribution to any sta- 
tistical analysis. Therefore, observations also do not support the 
idea of episodes of the Hall drift induced instabiUty in neutron 
stars. Note that recent calculations of the non-linear field decay 
by Pons & Geppert (2007) confirm the conclusion by HoUerbach 
& Riidiger (2002) that, even if the Hall-driven instability exists 



in neutron stars, it has no appreciable influence on the field de- 
cay. 

4. Statement of the problem 

To calculate the rate of Joule heating in young neutron stars, we 
consider a very simplified model. We assume that the magnetic 
field of neutron stars is generated by a turbulent dynamo at the 
beginning of their evolution. The dynamo a-parameter is propor- 
tional to the density gradient and, hence, the dynamo is most ef- 
ficient in the surface layers of proto-neutron stars (see Bonanno 
et al. 2005). Hence, one can expect that there exist electric cur- 
rents in the outer crust just after the convective stage. The decay 
of currents is highly inhomogeneous since the electric resistivity 
depends on the density and varies by many orders of magnitude 
over the crust. The decay timescale ists ~ A-kctL"^ /c^, where tr 
is the conductivity and L is the lengthscale of the magnetic field, 
and currents deposited closer to the surface will decay faster be- 
cause the conductivity is lower. Therefore, the thermal balance 
of young neutron stars with age f < 10^ yrs can be influenced 
only by dissipation of currents located in the surface layers. 

The decay of the crustal magnetic field is governed by the 
induction equation without the convective term 



dB 
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(5) 



where R = l/ir is the electric resistivity tensor and a is the 
conductivity tensor. In a magnetized crust, components of the 
conductivity tensor are 



Co ,cr_L 



O"0 

1 + 



(6) 



where (J|| and (7± are the components parallel and perpendicular 
to the magnetic field, an is the so-called Hall component that 
determines the current perpendicular to both the magnetic and 
electric fields, cto = e'^rieT/m^, is the coductivity at B = 0, 
Ue and are the number density of electrons and their effec- 
tive mass, ojb is the gyrofrequency of electrons, and t is the 
relaxation time. If the magnetic field is sufficiently strong and 
uibt > 1, then the conductivity components across the mag- 
netic field (T± and an are suppressed, and only the conductivity 
along the magnetic field is essential. Components of the resistiv- 
ity tensor R= 1 /a are given by 



R± = l/o-Q , Rh = wst/cto = B/ eerie 



(7) 



where is the number density of electrons. In a strong magnetic 
field, the components parallel and perpendicular to the magnetic 
field are the same and are not suppressed by the magnetic field. 
The Hall component of resistivity Rh can be higher than if 
ujbt > 1 that, generally, can be fulfilled in some layers of the 
neutron star. 

We neglect the Hall term in our consideration and, strictly 
speaking, our results apply to neutron stars with not very strong 
magnetic fields. However, as was argued in previous studies 
(Urpin & Shalybkov 1997; HoUerbach & Riidiger 2004; Pons 
& Geppert 2007), even in a very strong magnetic field that can 
magnetize plasma, the Hall effect insignificantly changes the rate 
of dissipation if the magnetic field contains a poloidal compo- 
nent. The difference in the rate of dissipation between the model 
that takes into account the Hall effect and the hypothetical model 
that has the same initial magnetic configuration but evolves in 
accordance with Eq. (5) without the Hall effect does not ex- 
ceed 10-20%. This is much smaller than uncertainties in our 
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knowlegde of the initial magnetic configuration or the chemi- 
cal composition of the crust. Such accuracy is sufficient for our 
purposes, and we will mimic the field decay by Eq. (5). 

We restrict our consideration to a dipolar field that can be de- 
scribed by the vector potential A = (0, 0, S{r, t) sm.9/r), where 
r and 9 are the radius and polar angle, respectively. The function 
5'(r, t) obeys the equation 

5^ _ 2S; _ ATTao_dS_ 

with the boundary conditions dS/dr + S/R = at the stellar 
surface r = R and S* ^ in the deep crustal layers. The cp- 
component of the electric current maintaining the dipolar mag- 
netic field is 

We have for the rate of Joule heating q{r, 9,t) = j^/ao. For 
simplicity, nonsphericity is often neglected in cooling calcula- 
tions. Therefore, we will use the polar-averaged expression for 
the rate of Joule heating. It is also convenient to normalize the 
function S{r, t) to its initial value at the surface, S{R, 0), which 
in turn can be related to the initial field strength at the magnetic 
equator by S{R, O) = R^B^. Then, the polar-averaged rate 
of Joule heating is given by 

._ c^R^B^, fd^s 2sY 

^ 247r2r2<To \dr^ J ' ^ ' 

In young neutron stars with age t < 10^ yrs, the conductivity 
of the outer crust and the outer part of the inner crust is deter- 
mined by electron scattering on phonons. Scattering on impuri- 
ties is important only in deep layers of the iimer crust even if 
the impurity parameter is large (see, e.g., Jones 2004). We use 
the conductivity due to electron-phonon scattering calculated by 
Itoh et al. (1993). These conductivities depend on the tempera- 
ture. Calculations presented here are based on the standard cool- 
ing scenario, which corresponds to a star with standard neutrino 
emissivities. We consider the field decay using the thermal his- 
tory for the 1 .4 Af0 neutron star constructed with the equation of 
state of Friedman & Pandharipande (1981) which is representa- 
tive of intermediate equations of state. Our choice is compelled 
by the fact that models with the equation of state stiffer than 
that of Friedman & Pandharipande agree better with data on the 
magnetic evolution of pulsars (Urpin & Konenkov 1997). 

Calculating the rate of Joule heating, we neglect its back in- 
fluence on the thermal evolution of a neutron star. It is not the 
purpose of this paper to discuss in detail how the thermal and 
magnetic evolution are coupled. One can neglect this coupling if 
the rate of heating does not exceed the luminosity of a standard 
neutron star. Our approach allows us to obtain the lower limit 
on the field strength that can affect the thermal evolution due to 
Joule heating, and we compare this limit to observational data. 
Note, however, that even the temperature of magnetars which 
are believed to be heated by the decay of extremely strong fields 
> 10^^ G is only a factor ~ 2 — 3 higher than the temperature 
of a standard cooUng neutron star. Therefore, our calculations 
can provide an order of magnitude estimate for the rate of Joule 
heating in the case of such strong magnetic fields. This simpli- 
fied approach can be useful in understanding the main qualita- 
tive features of Joule heating in magnetars. A study of the cou- 
pled early magnetothermal evolution of very strongly magne- 
tized stars such as magnetars is a more complicated problem, 
and it will be considered in a forthcoming paper. 



If the field in the surface layers is very strong, such as 
B"^ /A-K > /X where n is the shear modulus of the crust, then the 
induction and momentum equations are coupled, and magnetic 
stresses can induce motion of the crustal material. Therefore, de- 
cay of the crustal magnetic field can change the Lorentz force 
and, generaUy, can lead to crustal motion. The shear modu- 
lus is /U ~ 0.1ni{Ze)'^/a where a = (3/47rn)^/^, Z and rij 
is the charge and number density of ions (Strohmayer et al. 
1991). The condition B'^ /4n > /i is approximately equivalent 

to p < lO^Bi^^ g/cm3 where B13 = B/10^^ G, and one should 
expect that Joule heating of the surface layers in high magnetic 
field pulsars is accompanied by motion in the crust. We will ne- 
glect, however, the contribution of this effect to heating. 



5. Numerical results 

In calculations, we assume that the magnetic field generated by 
a dynamo initially occupies the outer layer of the crust up to a 
density po thus, initiaUy, B ^ at p < po- Since the field is 
confined to the crust, the density po can be lower or equal to 
the density at the crust-core boundary, ^ 2 x 10^'* g/cm''. It is 
difficult to determine the inner boundary of the generation re- 
gion from the physical data since both turbulent dynamo coeffi- 
cients and stability properties of proto-neutron stars are subject 
to many uncertainties. Therefore, we perform calculations for a 
few different values of pp. Note that the magnetic evolution is 
in good agreement with pulsar data if po ~ 10^^ g/cm^ (Urpin 
& Konenkov 1997). This conclusion is in qualitative agreement 
with the dynamo models of neutron stars since the dynamo is 
efficient only in the surface layers where the density gradient 
is large. After the dynamo action stops in a proto-neutron star, 
the field diffuses into deeper layers and can be non- vanishing at 
p > Po. The initial distribution of the magnetic field at p < po 
is also uncertain but, fortunately, the main qualitative results are 
not crucially sensitive to this distribution. We consider two mod- 
els of the initial distribution shown in Fig. 1 . These models differ 
by the characteristic depth of currents with model 2 (shown in 
the dashed line) corresponding to a smaller depth than model 1 
(solid line). 

In Fig. 2, we plot the distribution of the rate of Joule heating, 
q, for different ages. The distribution turns out to be complicated 
with a deep minimum that originates because, for a dipole field, 
the current density changes sign at some depth. The depth of 
this minimum increases with time due to diffusion of the mag- 
netic field into deep layers. The origin of this minimum is caused 
by the general diffusive properties of the crustal magnetic field 
and can be understood from the simplest quaUtative considera- 
tion. Indeed, if the initial distribution of s is like those shown 
in Fig. 1 , then the magnetic field will diffuse with time into deep 
layers with p > po and, hence, the time derivative of s will be 
positive somewhere at p > po, at least. On the another hand, 
the value of s and the field strength at the surface will obviously 
decrease with time because of ohmic dissipation and,hence, we 
have ds/dt < at r = R. Since s is a continuous function of the 
radius r and ds/dt has different signs at the surface and some- 
where in deep layers, there always exists at least one point inside 
the neutron star crust where ds/dt = 0. Combining Eqs. (8) and 
(9), we find that = at that point as well. Since q oc j^, the 
rate of heating is also vanishing at the same depth. In Fig. 2, we 
do not show q at this depth in detail because the region where 
g « is very narrow. The first maximum corresponds to a rela- 
tively low density ^ 10^" — 10^^ g/cm^, depending on the model. 
Initially, this maximum is located at a very low density, but it 
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Fig. 1. The dependence of the initial distribution of s(0, r) on the 
depth from the surface z = i? — r for the models 1 (solid line) 
and 2 (dashed line) and for po = 10^^ g/cm^ (the corresponding 
depth zo is approximately 330 m. 



moves slowly inward with the time. After 10^ yrs, the maximum 
reaches a depth ~ 300 m corresponding to p ~ 10^^ g/cm^. At 
that time, the distribution of q is already rather flat in the sur- 
face layers. The second maximum of q is located substantially 
deeper, at densities ~ 3 x 10^^ — 3 x 10^'^ g/cm^, and it also 
moves inward with time. The position of the second maximum 
is less sensitive to the initial models because our models 1 and 
2 differ mainly in the distribution of currenls in the surface lay- 
ers. Perhaps such a comphcated distribution of q can result in 
a comphcated radial dependence of the temperature in strongly 
magnetized stars. The rate of heating can reach a very high value 
^ lO'^^' — lO^'^ erg/cm'^-s or even higher if the initial magnetic 
field is stronger than S > 3 x lO^'' G. The total rates of heating 
above and below the region where g ~ are of the same order 
of magnitude. 

In Fig. 3, we plot the evolution of the density where the 
rate of Joule heating is maximal for models 1 and 2. This density 
corresponds to the first maximum of q (see Fig. 2). The location 
of the layer where the maxumum energy is released is important 
for understanding the energy budget in neutron stars (see, e.g., 
Kaminker et al. 2006). During the initial ^ lO"* yrs, the depth 
of this layer is very small and corresponds to < 10^^ g/cm'\ 
At this stage, depends on the initial distribution of currents. 
Later on, att > 10^ yrs, this depth becomes almost independent 
of the initial distribution. Up to t ~ 10* yrs, the maximum of 
Joule heating is located in the layers with a density lower than 
4 X 10^^ g/cm^. Likely, heat released in the layers with such 
a low density will be transfered to the surface and emitted by 
photons rather than diffusing to the core and being emitted by 
neutrinos. 

In Fig. 4, we show the evolution of the total rate of Joule 
heating integrated over the crust volume, Q = J qdV. For com- 
parison, we also plot in this figure the photon luminosity of a 
neutron star with standard cooling. If the initial magnetic field 
is strong and the corresponding heating rate is above the dotted 



igp 

Fig. 2. The dependence g on r for different ages and models. The 
sohd and dashed fines corresponds to the models shown in the 
sofid and dashed lines in Fig. 1. Curves 1, 2, and 3 show the 
distribution for t = 10^, 10^, and 10* yrs, respectively. 



fine, our results can provide only an order of magnitude esti- 
mate because Joule heating will influence the neutron star ther- 
mal evolution. The rate of Joule heating is calculated for three 
values of po = 10^^, lO^'', and 10^^ g/cnr^. For all these ini- 
tial distributions of the magnetic field, the rate of Joule heating 
can reach rather high values ~ 10"^^ — 10"^^ erg/s during ~ 10^ 
yrs if the initial field strength is 3 x 10^^ G. This heating rate is 
higher than the luminosity of a standard neutron star without any 
additional heating mechanisms. At a given initial field strength, 
the heating rate is initially higher for smaller po, but it decreases 
faster with the time. For example, if the magnetic field initiaUy 
occupies the surface layer with a density p < po = 10^^ g/cm^, 
then the heating rate at t ^ 100 yrs is approximately a factor 
of 10 higher than in the case of po = lO^"' g/cnr^. However, the 
heating rate for po = 10^^ g/cm^ decreases very rapidly and be- 
comes smaller already after ~ 2 x 10^ yrs. GeneraUy, the rate of 
Joule heating can be high in young neutron stars. Even in pulsars 
with initial magnetic fields ^ 3 x lO^'', Joule heating exceeds 
the photon luminosity during ~ 10'' yrs if po ~ 10^'^ g/cm^. 
Therefore, young pulsars with such magnetic fields should be 
hotter than it is predicted by the standard cooHng scenario. High 
magnetic field pulsars such as PSR J1847-0130 (B « 9 x 10^^ 
G) or PSR J1718-3718 (B « 7 x lO^^) can be in such a "heated" 
state even longer (~ 10^ yrs). However, after that time, their 
thermal evolution will follow the standard scenario. Our calcu- 
lations show that even the rate of heating required for magne- 
tars, Q ^ 10'*''' — 3 x 10''^ erg/s, is reachable during the initial 
evolution if the initial magnetic field at the surface is ~ 10^* 
G. Note that the rate of Joule heating depends not only on the 
field strength but also on the depth of the region where the mag- 
netic field was generated. This dependence is rather sensitive 
and, generally, it is possible that the heating rate differs by more 
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igt 

Fig. 3. The depth of the first maximum of g as a function of time 
for the models shown in Fig. 1. 



than an order of magnitude in pulsars of the same age and hav- 
ing approximately the same surface magnetic fields (for exam- 
ple, compare the soUd and dash-and-dotted lines in Fig. 4). Such 
different heating rates can lead to a substantial difference in the 
surface temperature in high magnetic field pulsars. 



Igt 

Fig. 4. The dependence of Q on time for different models. The 
magnetic field is normahzed to 3 x 10^'^ G. The sohd and dashed 
lines correspond to the initial models 1 and 2 from Fid. 1, re- 
spectively, with po = 10^'^ g/cm'^. The dashed-and-dotted and 
dotted-and-dashed hues show Q for the model 1 with po = 10^^ 
and 10^^ g/cm^, respectively. The dotted line shows the time de- 
pendence of the luminisity of a neutron star with no additional 
heating mechanism. 



6. Discussion 

We have considered dissipation of a magnetic field initially an- 
chored to a neutron star crust. The magnetic field of neutron stars 
can be generated by a turbulent dynamo during a short hydro- 
dynamically unstable phase at the beginning of their evolution. 
The mean-field dynamo is most effecient in the surface layers 
where the density gradient is large and, hence, the dynamo a- 
parameter is maximal. Therefore, electric currents maintaining 
the magnetic configuration are likely located in the surface lay- 
ers. Dissipation of electric currents is a highly inhomogeneous 
process because the density and, hence, conductivity varies by 
many orders of magnitude through the crust. The rate of dissipa- 
tion decreases with increasing density and, hence, the field decay 
is slower in deep crustal layers. Dissipation of currents located 
in the surface layers proceeds on a relatively short time-scale and 
is accompanied by massive heat release. Joule heating is so ef- 
ficient that it can even influence the thermal evolution of young 
neutron stars if the magnetic field is sufficiently strong. 

The duration of the stage when Joule heating can affect the 
thermal evolution is determined by the strength of the magnetic 
field and the depth where this field is anchored. Both these quan- 
tities can vary over a wide range because they depend, for exam- 
ple, on poorly known details of the generation process and post- 
collapse hypercritical accretion that can essentially submerge 
the generated magnetic field (Chevalier 1989). At a given field 
strength, the rate of Joule heating is higher in a star where the 
field was initially anchored in the layer with a lower density. 



However, the time interval when this heating is important for the 
thermal evolution is shorter in such stars. On the contrary, if the 
field is anchored in a deeper layer, Joule heating is less effecient 
but it can contribute significantly in the thermal evolution over 
longer time. 

Since the rate of Joule heating depends on two parameters it 
is generally possible that neutron stars with the same age and sur- 
face magnetic field can have substantially different heating rates 
and, hence, surface temperatures. This allows us to understand 
why the surface temperature of the high magnetic field pulsar 
PSR J17 18-37 18 is substantially lower than that of magnetars 
despite its magnetic field being approximately that of magnetars 
(7.4 X 10^^). The magnetic field of high magnetic field pulsars 
was likely anchored initially in layers with a lower density, and 
the stage of efficient Joule heating when their luminosity was 
comparable to that of magnetars is completed. 

Our calculations show that additional heating provided by 
dissipation of the crustal magnetic field can even be comparable 
to the luminosity of magnetars (~ 10'^^ — 10'^^ erg/s) if the mag- 
netic field at the stellar surface is as high as ~ (1 — 3) x 10^^ 
G. Generally, such a strong magnetic field can be generated in 
rapidly rotating proto-neutron stars by the turbulent a-dynamo 
(Bonanno et al. 2006). Note that our calculations provide only 
the lower limit on the rate of Joule heating since we neglected 
the back influence of this heating on the crustal temperature. 
Additional heating will increase the temperature and decrease 
the conductivity and, hence, it will result in a higher dissipation 
rate if one consistently takes into account the coupling between 
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Table 1. 

The surface temperature and magnetic field of young radiopulsars. 
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Fig. 5. The surface temperature Tg versus the magnetic field for 
22 young neutron stars. 



the thermal evolution and field decay. For strongly magnetized 
stars with a surface field ^ (1 — 3) x 10^^ G, the epoch when 
Joule heating is of the order of the magnetar luminosity can last 
as long as ~ 10* — 10^ yrs if the magnetic field was initially 
anchored in the layers with pQ ~ 10^^ g/cm^. This time is com- 
parable to the estimated life-time of magnetars. If the initial field 
is weaker ^ 3 x 10^'^ G, but is confined to the same layers, then 
the rate of Joule heating exceeds the standard luminosity of a 
cooUng neutron star for a shorter time, 10^ yrs. During this 
time, the rate of heating is approximatelly one order of magni- 
tude greater than the luminosity of a standard neutron star with 
the same age. Even in neutron stars born with a relatively weak 
magnetic field ~ 3 x 10^^ G, Joule heating can influence the 
early thermal elolution but for a very short time and only if the 
magnetic field was initially confined to low density layers. 

From our results, we conclude that many young neutron stars 
can ehxibit departures from the standard cooling scenario if the 
magnetic field is of crustal origin and if it is sufficiently strong. 
How large these departures are and how long Joule heating can 
manifest is determined by the initial strength of the magnetic 
field and its location. At the moment, observations only provide 
information on the surface temperature of neutron stars with an 
age greater than ~ 1 kyr. For stars with such age, our model pre- 
dicts that departures from standard cooling can be expected only 
if the magnetic field is sufficiently strong, B > (1 — 3) x 10^^ G. 
Therefore, recently discovered high magnetic field pulsars such 
as PSR J1847-0130 (9.4 x lO^^ G), PSR J1718-3718 (7.4 x lO^^ 
G), PSR J1814-1744 (5.5 x lO^^ G), PSR Jl 1 19-6127 (4.4x lO^^ 
G), and PSR B0154H-61 (2.1 x lO^^ G) all can have a higher sur- 
face temperature than is predicted by standard cooUng models 
without additional heating. 

In Fig. 5, we plot the available data on the surface tempera- 
ture Ts and magnetic fields B of 22 young neutron stars. These 
data include 9 magnetars listed in Table 1 by Pons et al. (2007) 



Source 


kT 


B 


Ref. 




(keV) 


(10^2 G) 




RX J0822-4320 


0.14-0.16 


3.4 


Zavlinetal.(1999) 








Pavlov et al.(1999) 


IE 1207.4-5209 


0.12-0.16 


0.35 


Zavlin et al.(2004) 








Gotthelf etal.(2007) 


Vela 


0.06-0.07 


3.4 


Page et al.(2004) 


PSR B 1706-44 


0.07 


3.0 


Page et al.(2004) 


PSRJ0538-I-2817 


0.075 


0.8 


Levandowski et al. 








(2004) 


Geminga 


0.043 


1.6 


De Luca et al.(2005) 


PSR 1055-52 


0.064 


1.15 


De Luca et al.(2005) 


PSRB0656+14 


0.077 


4.7 


De Luca et al.(2005) 


PSRB2334+61 


0.056 


10 


McGowan et al. 








(2006) 


PSRJ1718-3718 


0.145 


74 


Kaspi at al.(2005) 


PSR J1119-6127 


0.206 


41 


Gonzalez et al.(2005) 


CXOU J1819-1458 


0.120 


50 


Reynolds et al.(2006) 


PSR J 1357-6429 


0.085 


8 


Zavlin (2007) 



which are likely young objects according to current theory. We 
also plot the data on B and Tg for 13 radiopulsars with a spin- 
down age shorter than 1 Myr If the magnetic fields of these pul- 
sars decay, their true age should be substantially less than the 
spin-down age and, Ukely, such pulsars are sufficiently young 
to be treated in our analysis. These 13 radiopulsars are listed in 
Table 1. 

In accordance with our calculations, the thermal history of 
young pulsars with an age < 10^ yrs can be influenced by 
Joule heating if the magnetic field at birth is stronger than 
(1 — 3) X lO^'"* G. Indeed, it seems that the surface temperature of 
high magnetic field pulsars with B > 10^^ G shows some trend 
to be higher than that of low magnetic field pulsars, in qualita- 
tive agreement with the prediction of our model. This difference 
cannot be accounted for by the influence of the magnetic field 
on the thermal conductivity, because the geometry and strength 
of the magnetic field are unimportant for the surface averaged 
thermal structure of neutron stars if the field strength is moder- 
ate, B < 10" G (Potekhin et al. 2005). For such B, the mag- 
netic field can affect the thermal evolution of neutron stars only 
via Joule heating, and the surface temperature should be propor- 
tional to approximately B^/^. The data in Fig. 5 obviously are 
not sufficient to prove this dependence but, at least, they do not 
contradict it. More detections of the surface temperature in high 
magnetic field pulsars are needed, and they could provide strong 
evidence of the crustal origin of the magnetic field in neutron 
star 

It is important to consider coupled magnetothermal evolu- 
tion by solving the heat and induction equations because char- 
acteristic timescales of the thermal and crustal magnetic evolu- 
tion are approximately comparable in young neutron stars with 
t < 10^ yrs. However, this problem, taking into account the 
magnetically-induced anisotropy of the thermal conductivity and 
electric resistivity tensors,is far from being solved and the study 
of particular problems of the magnetothermal evolution can help 
in understanding when Joule heating can affect neutron star cool- 
ing. From the result of the present study and conclusions ob- 
tained by Miralles et al. (1998), Joule heating can be important 
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only in very young (t < 10'' - 10^ yrs) or in old (t > 10 Myr) ZavUn, v., Truemper, J., & Pavlov, G. 1999, ApJ, 525, 959 
neutron stars provided the magnetic field is sufficiently strong 
at birth. In stars with a moderate age, Joule heating seems to be 
inefficient. 
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